Centrality dependence of charged hadron transverse momentum spectra 
in Au+Au collisions from ^/s^ = 62.4 to 200 GeV 
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We have measured transverse momentum distributions of charged hadrons produced in Au+Au 
collisions at y / s^" = 62.4 GeV. The spectra are presented for transverse momenta 0.25 < pr < 
4.5 GeV/c, in a pseudo-rapidity range of 0.2 < -q < 1.4. The nuclear modification factor Raa is 
calculated relative to p+p data at the same collision energy as a function of collision centrality. 
For pr > 2 GeV/c, Raa is found to be significantly larger than in Au+Au collisions at ^/s NN = 
130 and 200 GeV. In contrast, we find that the evolution of the invariant yields per participant 
pair from peripheral to central collisions is approximately energy-independent over this range of 
collision energies. This observation challenges models of high pr hadron suppression in terms of 
parton energy loss. 

PACS numbers: 25.75.-q,25.75.Dw,25.75.Gz 



The yield of charged hadrons produced in collisions of 
gold nuclei at an energy of y^s"™ — 62.4 GeV has been 
measured with the PHOBOS detector at the Relativistic 
Heavy Ion Collider (RHIC) at Brookhaven National 
Laboratory. The data are presented as a function of 
transverse momentum (pt) and collision centrality. The 
goal of these measurements is to study the modification 
of particle production in the presence of the produced 
medium by comparing to nucleon-nucleon collisions at 
the same energy. 

This measurement was motivated by results from 
Au+Au collisions at y/s^ = 130 and 200 GeV. Hadron 
production at these energies was found to be strongly 
suppressed relative to expectations based on an indepen- 
dent superposition of nucleon-nucleon collisions at p? of 
2-10 GeV/c 0, H H 0. The modification of high-p T 
hadron yields has commonly been investigated using the 
nuclear modification factor Raa, defined as 

R AA - a pp el d 2N AA/dp T dri 
(Ncoii) d 2 a pp /dp T dri 

A value of Raa — 1 corresponds to scaling of particle 



production with the average number of binary nucleon- 
nucleon collisions, (N co u), within a heavy-ion collision. 

For the production of charged hadrons in central 
Au+Au collisions at v / s^" = 200 GeV, values of Raa ~ 
0.2 are observed at p T = 4 GeV/c H 0, 0. Such a 
suppression had been predicted to occur as a consequence 
of the energy loss of high-p*r partons in the dense medium 
formed in Au+Au collisions This hypothesis is 

consistent with the observed absence of this effect in 
deuteron-gold collisions at the same collision energy 




While studies of the system-size dependence of the 
high-pT suppression exist, the dependence on the colli- 
sion energy is poorly known. Data on the production of 
neutral pions in Pb+Pb collisions at y/s NN — 17.2 GeV 
suffer from large uncertainties in the parametrizations of 
the p+p reference data |Rlllllll2| . The results presented 
here for Au+Au collisions at y/s NN — 62.4 GeV, in 
combination with the results from 130 and 200 GeV 
collisions, allow the first detailed examination of the 
connection between high px suppression and collision 
energy, charged-particle density and collision geometry. 
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FIG. 1: Invariant yields for charged hadrons from Au+Au 
collisions at y/s NN = 62.4 GeV, in the pseudo-rapidity interval 
0.2 < r) < 1.4 as a function of pr for 6 centrality bins. 
For clarity, consecutive bins are scaled by factors of 2. The 
systematic uncertainties are smaller than the symbol size. 



The data were collected using the PHOBOS two- 
arm magnetic spectrometer. Details of the experimental 
setup can be found in |l3j . The primary event trigger 
used the time difference between signals in two sets of 
10 Cerenkov counters, located at 4.4 < \r/\ < 4.9, to 
select collisions along the beam-axis that were close to 
the nominal vertex position. 

For the analysis presented here, events were divided 
into six centrality classes, based on the observed signal 
in two sets of 16 scintillator counters covering pseudo- 
rapidities 3.0 < \rj\ < 4.5. The results of the Glauber 
calculation implemented in the HIJING model 0, 0] 
were used to estimate the average number of participat- 
ing nucleons, (N part }, and binary collisions, (N co u), for 
each centrality class. For the Glauber calculation, as 
well as for the determination of Raa at 62.4 GeV, we 
mb. The values and systematic 
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uncertainties for (N part ) and (N co u) at both 62.4 and 
200 GeV are listed in Table I. 

The event selection and track reconstruction procedure 
for this analysis closely follows the procedure for the 
previously published analysis at y / s^" = 200 GeV Q. 
Events with a primary vertex position within ±10 cm of 
the nominal vertex position were selected. Only particles 
traversing a full spectrometer arm were included in the 
analysis, resulting in a low transverse momentum cutoff 



TABLE I: Details of the centrality classes used in this 
analysis. Bins are expressed in terms of percentage of the 
total inelastic Au+Au cross-section. 



Centrality 


/N 62A \ 

\ J v part / 


(N™it) 


(N 200 <) 

\ J ' part 1 




45-50% 


61 ±7 


76 ± 12 


65 ±4 


107 ±16 


35-45% 


86 ±9 


120 ±16 


93 ±5 


175 ± 25 


25-35% 


130 ± 10 


215 ±22 


138 ±6 


300 ± 39 


15-25% 


189 ±9 


370 ± 24 


200 ±8 


500 ± 60 


6-15% 


266 ±9 


590 ± 21 


276 ±9 


780 ± 86 


0-6% 


335 ±11 


820 ± 25 


344 ± 11 


1050 ± 105 
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FIG. 2: Nuclear modification factor Raa for central A+A 
events, as a function of collision energy for pr = 1, 2, 3 
and 4 GeV/c. Filled symbols show data for charged hadrons, 
open symbols show tv° data. For the WA98 7r° data, Raa 
is shown using three different parameterizations for the p+p 
reference spectrum [Tot ITU H2I : for pr — 4 GeV/c, the arrow 
indicates the lower limit on Raa- The error bars show the 
combined systematic and statistical uncertainty obtained by 
interpolating the pr dependence of Raa- 



at p T w 0.2 GeV/c. 

The transverse momentum distribution for each cen- 
trality bin was corrected separately for the geometri- 
cal acceptance of the detector, the inefficiency of the 
tracking algorithm, secondary and incorrectly recon- 
structed particles, and the distortion due to binning 
and momentum resolution. The relative importance of 
most of the corrections and the associated estimated 
systematic uncertainties are similar to those reported in 
the previous analysis 4]. The exception is the correction 
for momentum resolution and binning effects at high px, 
where changes to the track fitting procedure improved 
the resolution and reduced the corresponding uncertainty 
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FIG. 3: Ratio of pr distributions from Au+Au collisions to various reference distributions at yJs NN — 62.4 GeV (filled symbols) 
and 200 GeV (open symbols). Data are shown in six bins of centrality, ranging from {N par t) = 61 to 335 for 62.4 GeV collisions. 
The top row shows the nuclear modification factor Raa, i.e. the ratio relative to proton-(anti)proton collisions scaled by (N co ii)- 
The middle row shows R A F A art , which uses proton-(anti)proton spectra scaled by (N part /2). The bottom row shows R P Q art , 
using a fit to central data scaled by {N par t/2) as a reference. The dashed line in the middle and bottom rows indicates the 
expectation for (N co u) scaling at ^/s NN = 62.4 GeV relative to the reference distribution. Systematic uncertainties for all plots 
are shown by brackets (90% C.L.) 



by a factor of two. 

In Fig. ^ we present the invariant yield of charged 
hadrons as a function of transverse momentum, obtained 
by averaging the yields of positive and negative hadrons. 
Data are shown for six centrality bins and are averaged 
over a pseudo-rapidity interval 0.2 < rj < 1.4. 

Focusing first on central events, Fig. [21 shows the 
nuclear modification factor Raa for the 6% most central 
Au+Au collisions at y/s NN = 62.4 GeV for four different 
values of pt ranging from 1 to 4 GeV/c. Results from 
p+p collisions at the same energy were used in the 
calculation of Raa- The new results are compared with 
results for charged hadrons at y 7 ?^" = 130 and 200 GeV. 
For pt — 2 GeV/c and above, we observe a smooth 
decrease of Raa from 62.4 to 200 GeV in central Au+Au 
collisions. 

Also shown in Fig. [21 is Raa for ti -0 production in 
Pb+Pb collisions at 17.2 GeV from WA98, obtained 



using three p+p reference parametrizations 0, 0> EH j 
as well as tt° data for 130 GeV and 200 GeV from 
PHENIX [H|l3. Although the uncertainties in the 17.2 
GeV reference distribution are large, the data indicate 
that Raa for n° production at pt > 3 GeV/c drops 
from Raa > 1 at ^/s^ = 17.2 GeV to Raa < 0.2 
at y/s NN = 200 GeV. It is interesting to note that the 
change in the observed charged hadron pseudo-rapidity 
density over the same energy range is only a factor of two 
[19) . The data also show that Raa for neutral pions is 
consistently lower than for charged hadrons at the same 
collision energy, which is important when comparing data 
from different experiments. 

The centrality evolution of Raa for the 62.4 GeV data 
can be studied in detail in the top row of Fig. [3 For com- 
parison, we include the results from Au+Au collisions 
at 200 GeV 4], using the same centrality binning. We 
observe that Raa reaches a maximum of approximately 
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1.6 for peripheral collisions at 62.4 GeV, whereas the 
maximum value observed at 200 GeV is close to 1. For 
both datasets, a significant and systematic decrease in 
Raa is observed in progressing from peripheral to central 
events. 

It has been previously noted that the observed strong 
centrality dependence of Raa at ^/s^ = 200 GeV 
corresponds to a relatively small change in the yield 
per participating nucleon Q and that over the same 
centrality range, the total yield of charged particles per 
participating nucleon is constant within the experimental 
uncertainties [2(|. These observations lead us to define 

R N part _ gpp e ' d 2 N A A/dp T d7] 
AA (N part /2) d 2 a pp /dp T d V ' 1 > 

in analogy to Equation 1, where now we scale the 
reference spectrum by N part /2, rather than N co u. The 

centrality, px and energy dependence of R AA art is shown 
in the middle row of Fig.[3]for collisions at 62.4 GeV and 
200 GeV. Over the full range of px studied here, the yield 
per participant pair shows a variation of only rs 25% from 
peripheral to central collisions for both collision energies. 

This centrality independence is further illustrated in 
the bottom row of Fig. [3J The quantity Rp^ rt , defined 
as 

R N part _ ( N parf°) d 2 N A A/dp T dr] 
PC (N part ) d 2 N° A - A 6% /dp T dv [ ) 

is shown as a function of px for the six centrality bins. 
R p P Q Tt measures the change in yield per participant pair, 
relative to a fit to the central data. Normalizing to 
the central, rather than the peripheral, data has the 
advantage that the results are more easily compared 
among different experiments with different ranges in 
centrality, while providing the same information in com- 
parisons with theoretical calculations. Data are shown 



for collisions at 62.4 GeV and 200 GeV. This plot again 
shows the small variation of the yield per participant 
pair from peripheral to central collisions. Furthermore, 
it demonstrates that the modification of the yield from 
peripheral to central collisions is the same for both 
energies, over the full px range, within experimental 
uncertainties of less than 10%. This striking agreement 
can be compared with the much larger variation of Raa 
as a function of energy, centrality and px- 

This measurement of the centrality dependence of 
particle yields at yfs~^ = 62 .4 GeV severely constrains 
the parameter space for any model of high-py hadron 
suppression in terms of parton energy loss. It remains 
to be seen whether the agreement in the centrality 
evolution of the yields per participant pair, which ranges 
over a factor of three in collision energy, is a natural 
consequence of such models. It is important to note 
that the energy-independence of the centrality evolu- 
tion is a characteristic feature , no t only of total and 
differential particle yields [2(], l2l| . but also of multi- 
particle correlation measurements such as Bose-Einstein 
correlations p2| . The apparent dominance of the initial- 
state geometry, even for observables closely related to 
the dynamical evolution of heavy-ion collisions, is one of 
the key features of these interactions that remains to be 
understood. 
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